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Abstract 

Black branes in AdSs appear in a four parameter family labeled by their velocity 
and temperature. Promoting these parameters to Goldstone modes or collective 
coordinate fields - arbitrary functions of the coordinates on the boundary of AdSs 
- we use Einstein's equations together with regularity requirements and boundary 
conditions to determine their dynamics. The resultant equations turn out to be 
those of boundary fluid dynamics, with specific values for fluid parameters. Our 
analysis is perturbative in the boundary derivative expansion but is valid for arbitrary 
amplitudes. Our work may be regarded as a derivation of the nonlinear equations of 
boundary fluid dynamics from gravity. As a concrete application we find an explicit 
expression for the expansion of this fluid stress tensor including terms up to second 
order in the derivative expansion. 
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1 Introduction 

The AdS / CFT correspondence provides an important laboratory to explore both gravita- 
tional physics as well as strongly coupled dynamics in a class of quantum field theories. 
Using this correspondence it is possible to test general lore about quantum field theory in 
a non perturbative setting and so learn general lessons about strongly coupled dynamics. 
Conversely, it is also possible to use the AdS / CFT duality to convert strongly held convic- 
tions about the behaviour of quantum field theories into general lessons about gravitational 
and stringy dynamics. 

In this paper we use the AdS/CFT correspondence to study the effective description 
of strongly coupled conformal field theories at long wavelengths. On physical grounds it 
is reasonable that any interacting quantum field theory equilibrates locally at high enough 
energy densities, and so admits an effective description in terms of fluid dynamics. The 
variables of such a description are the local densities of all conserved charges together with 
the local fluid velocities. The equations of fluid dynamics are simply the equations of local 
conservation of the corresponding charge currents, supplemented by constitutive relations 
that express these currents as functions of fluid mechanical variables. As fluid dynamics 
is a long wavelength effective theory, these constitutive relations are usually specifled in a 
derivative expansion. At any given order, thermodynamics plus symmetries determine the 
form of this expansion up to a flnite number of undetermined coefficients. These coefficients 
may then be obtained either from measurements or from microscopic computations. 

The best understood examples of the AdS/CFT correspondence relate the strongly 
coupled dynamics of certain conformal field theories to the dynamics of gravitational sys- 
tems in AdS spaces. In this paper we will demonstrate that Einstein's equations with 
a negative cosmological constant, supplemented with appropriate regularity restrictions 
and boundary conditions, reduce to the nonlinear equations of fiuid dynamics in an ap- 
propriate regime of parameters. We provide a systematic framework to construct this 
universal nonlinear fiuid dynamics, order by order in a boundary derivative expansion. 
Our work builds on earlier derivations of linearized fiuid dynamics from linearized gravity 
by Policastro, Son and Starinets and on earlier examples of the duality between non- 
linear fiuid dynamics and gravity by Janik, some of the current authors and collaborators 

I, H, H (c/, [|, ^ H, I, [To], [ill for related work and extensions and fl^, |l^ for some recent 



work). There is a large literature in deriving linearized hydrodynamics from AdS/CFT, 



see [y, |T5l||, |T8|, Ig, H, ||, |2|, |5|, H, H, gg, P for developments in 



this area and |32[ for a review and comprehensive set of references. 



Our results, together with those of earlier papers referred to above, may be interpreted 
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from several points of view. First, one may view them as a confirmation that fluid dy- 
namics is the correct long wavelength effective description of strongly coupled field theory 
dynamics. Second, one could assume the correctness of the fluid description and view our 
results as providing information on the allowed singularities of 'legal' solutions of gravity. 
Finally, our work may be used to extract the values of all coefficients of the various terms 
in the expansion of the stress tensor in the fluid dynamical derivative expansion, for the 
fluid dual to gravity on AdSs. The universal behaviour of the shear viscosity - a coefficient 
of a term in the expansion of the stress tensor to first order in field theory derivatives - in 
fluids dual to gravity ||2^ has already attracted attention and has impacted experimental 



analysis of RHIC data ||33| , 0, In this paper we work out the universal values of all 
coefficients of (nonlinear) two derivative terms stress tensor of the distinguished conformal 
fluid dual to gravity on AdSs. 

Consider any two derivative theory of five dimensional gravity interacting with other 
fields, that has AdSs as a solution. Examples of such theories include IIB supergravity 
on AdSs where M is any compact five dimensional Einstein manifold with positive 
cosmological constant; for example M = S^, T^'^ and Y^'"^ for all p, q. The solution space 
of such systems has a universal sub-sector; the solutions of pure gravity with a negative 
cosmological constant.[| We will focus on this universal sub-sector in a particular long 
wavelength limit. Specifically, we study all solutions that tubewiseQ approximate black 
branes in AdSs, whose temperature and boost velocities vary as a function of boundary 
coordinates x'^ on a length scale that is large compared to the inverse temperature of the 
brane. We investigate all such solutions order by order in a perturbative expansion; the 
perturbation parameter is the length scale of boundary variation divided by the thermal 
length scale. Within the domain of validity of our perturbative procedure (and subject 
to a technical assumption), we establish the existence of a one to one map between these 
gravitational solutions and the solutions of the equations of a distinguished system of 
boundary conformal fluid dynamics. Implementing our perturbative procedure to second 
order, we explicitly construct the fluid dynamical stress tensor of this distinguished fluid 
to second order in the derivative expansion. 

Roughly speaking, our construction may be regarded as the 'Chiral Lagrangian' for 
brane horizons. Recall that the isometry group of AdSs is 5*0(4, 2). The Poincare algebra 
plus dilatations form a distinguished subalgebra of this group; one that acts mildly on 
the boundary. The rotations SO (3) and translations R^'^ that belong to this subalgebra 
annihilate the static black brane solution in AdSs. However the remaining symmetry gen- 



^ Recall that the Einstein frame Lagrangian contains no interaction terms that are linear in the non 
gravitational fluctuations. 

^We will work in AdS spacetimes where the radial coordinate r S (0, oo) and will refer to the remaining 
coordinates = {v,Xi) € R^'"^ as field theory or boundary coordinates. The tubes referred to in the text 
cover a small patch in field theory directions, but include all values of r well separated from the black 
brane singularity at r = 0; typically r > where is the scale set by the putative horizon. 
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erators - dilatations and boosts - act nontrivially on this brane, generating a 4 parameter 
set of brane solutions. These four parameters are simply the temperature and the velocity 
of the brane. Our construction effectively promotes these parameters to 'Goldstone fields' 
(or perhaps more accurately collective coordinate fields) and determines the effective dy- 
namics of these collective coordinate fields, order by order in the derivative expansion, but 
making no assumption about amplitudes. Of course the collective coordinates method has 
a distinguished tradition in theoretic physics; see for instance the derivation of the Nambu- 
Goto action in Our paper, which applies these methods to black brane horizons, is 



strongly reminiscent of the membrane paradigm of black hole physics, and may perhaps be 
regarded as the precise version of this paradigm in its natural setting, i.e., AdS spacetime. 

Seen from inverse point of view, our construction may be regarded as a map from 
solutions of the relativistic fluid dynamics equations on R^'^ to the space of long wavelength, 
locally black brane, solutions of gravity in AdSs. That is, we present a systematic procedure 
to explicitly construct a metric dual to any solution of the equations of the distinguished 
fluid dynamics alluded to above. This metric solves the Einstein's equations to a given 
order in the derivative expansion (one higher than the order to which the equations of fluid 
dynamics were formulated and solved), asymptotes to AdSs with a boundary stress tensor 
equal to the fluid dynamical stress tensor, and is regular away from the usual singularity 
of black branes (chosen by convention to be at r = 0). 

As an important physical input into our procedure, we follow |^, ^ to demand 
that all the solutions we study are regular away from the r = curvature singularity of 
black branes, and in particular at the the location of the horizon of the black brane tubes 
out of which our solution is constructed. We present our construction in the analogue 
of Eddington-Finklestein coordinates which extend all the way to the future curvature 
singularity. Although we have not yet performed a careful global analysis of our solutions, it 
seems rather clear that they each possess a regular event horizon that shields the boundary 
from this curvature singularity. 

This paper is organized as follows. We begin in § |2| with the basic outline of the 
computation expanding on the ideas presented above. In § ^ we outline in detail the logic 
and strategy of our perturbative procedure. We then proceed in § | to implement our 
perturbative procedure to first order in the derivative expansion. In § ^ we extend our 
computation to second order in the same expansion. In § ^ we demonstrate the Weyl 
invariance of the fluid dynamical stress tensor we obtain, and further use this stress tensor 
to compute corrections to the dispersion relation for sound and shear waves in this fluid. 
In § 1^ we end with a discussion of our results and of future directions. 

Note added: After we had completed writing this paper we learnt of related work soon 



to appear |33]. The authors of this paper utilize Weyl invariance to constrain the form of 



the second order fluid dynamical stress tensor up to 5 undetermined coefficients. They then 
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use information from linearized gravitational quasinormal mode calculations together with 
an earlier computation of Janik and collaborators to determine 3 of these five coefficients. 
As far as we have been able to tell, their results are consistent with the full second order 
stress tensor (and prediction for quasinormal mode frequencies) presented herein. This is a 
nontrivial check of our results. We thank the authors of ||3^ for sharing their results with 
us prior to publication. 

Note added in v2: The preprint appeared simultaneously in an arXiv listing with vl 



of this paper. In § ^]3| of this updated version of our preprint we have presented a detailed 
comparison of our results with those of p^; where they overlap we find perfect agreement. 



2 Fluid dynamics from gravity 

We begin with a description of the procedure we use to construct a map from solutions 
of fiuid dynamics to solutions of gravity. We then summarize the results obtained by 
implementing this procedure to second order in the derivative expansion. 

Consider a theory of pure gravity with a negative cosmological constant. With a par- 
ticular choice of units {RAds = 1) Einstein's equations are given by0 

Emn = Rmn — -9mnR — 6 9mn = {2V) 
=^ Run + 4 guN = 0, R = -20. 

Of course the equations ( |2.1| ) admit AdSs solutions. Another class of solutions to these 
equations is given by the 'boosted black branes'0 

ds^ = -2 Uf, dx^'dr - f{hr) dx^'dx" + P^^ dx'^dx" , (2.2) 



with 



fir) = l-\ 



J.4 



U" 



U 



1 



■^We use upper case Latin indices {M, A^, • • • } to denote bulk directions, while lower case Greek indices 
{/^, v^- ■ ■} refer to field theory or boundary directions. Finally, we use lower case Latin indices {«, j, • • • } 
to denote the spatial directions in the boundary. 

*The indices in the boundary are raised and lowered with the Minkowski metric i.e., = rj^^ u'^ . 
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where the temperature T = ^ and velocities j3i are all constants with = (3j (3\ and 

p/^^ = + 7]^"" (2.4) 



is the projector onto spatial directions. The metrics ( p.2|) describe the uniform black brane 
written in ingoing Eddington-Finkelstein coordinates, at temperature T, moving at velocity 

Now consider the metric ( p.2|) with the constant parameter b and the velocities /3j 
replaced by slowly varying functions h{x^) ^ (5i{x^) of the boundary coordinates. 

ds^ = -2 rfx^ dr-r'^f r) rfx'^ rfx"^ + P^^(x") rfx^' c/x'^ . (2.6) 

Generically, such a metric (we will denote it by g^^\h{x^) , (3i{x^)) is not a solution to 
Einstein's equations. Nevertheless it has two attractive features. Firstly, away from r = 0, 
this deformed metric is everywhere non-singular. This pleasant feature is tied to our 
use of Eddington-FinkelsteinQ coordinates.]] Secondly, if all derivatives of the parameters 
h{x^) and Pi{x^) are small, g^'^^ is tubewise^ well approximated by a boosted black brane. 
Consequently, for slowly varying functions b{x^), (3i{x^), it might seem intuitively plausible 
that (|2^ ) is a good approximation to a true solution of Einstein's equations with a regular 
event horizon. The main result of our paper is that this intuition is correct, provided the 
functions h{x^) and Pi{x^) obey a set of equations of motion, which turn out simply to be 
the equations of boundary fluid dynamics. 

Einstein's equations, when evaluated on the metric g^^\ yield terms of first and second 
order in field theory [i.e., (xj, v) = x^) derivatives of the temperature and velocity fields.^ 
By performing a scaling of coordinates to set h to unity (in a local patch), it is possible to 
show that field theory derivatives of either \n.h{x^) or (3i{x^) always appear together with a 
factor of h. As a result, the contribution of n derivative terms to the Einstein's equations is 
suppressed (relative to terms with no derivatives) by a factor of (b/L)"' ~ 1/(TL)'". Here L 



As we have explained above, the 4 parameter set of metrics (2.2) may ah be obtained from 



ds^ =2dvdr--r^ f{r) dv^ + dx^ , (2.5) 

with f = ^ — ^ via a coordinate transform. The coordinate transformations in question are generated by 
a subalgebra of the isometry group of AdSs. 

^It is perhaps better to caU these gencraUzed Gaussian null coordinates as they are constructed with 
the aim of having the putative horizon located at the hypersurface r{x^) — r^- 

^A similar ansatz for a black branes in (for instance) Feffcrman- Graham coordinates i.e., Schwarzschild 
like coordinates respecting Poincare symmetry, is singular at rb= I. 

®As explained above, any given tube consists of all values of r well separated from ?■ = 0, but only a 
small region of the boundary coordinates x^. 

^As (J*-"-' is an exact solution to Einstein's equations when these fields are constants, terms with no 
derivatives are absent from this expansion. 
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is the length scale of variations of the temperature and velocity fields in the neighbourhood 
of a particular point, and T is the temperature at that point. Therefore, provided L T ^ 1, 
it is sensible to solve Einstein's equations perturbatively in the number of field theory 
derivatives -p^ 

In § we formulate the perturbation theory described in the previous paragraph, and 
explicitly implement this expansion to second order in 1/{LT). As we have mentioned 
above it turns out to be possible to find a gravity solution dual to a boundary velocity and 
temperature profile only when these fields obey the equation of motion 



d^Tf"" = 



(2.7) 



where the rescaleclQ stress tensor T^*^ (to second order in derivatives) is given by 
T^"" ={71 {rf" + 4m'^m^) - 2 {jxTf a^"" 

+ (7rT)2 [(In 2) + 2 T,^ + (2 - In 2) V- Ti^: + TZ + 



(2.8) 



where 



a 



rpflU 

-'■2a 
-'2d 



^2e 



1 

~ 3 



P^" a^'f^a. 



(2.9) 



Pn^ Vu'' - - Pf"" Vu^ 
3 

p,a p ^Q^^^^^^ _ 1 pa/3 p 



U 



Our conventions are 60123 = — e*^^^^ = 1 and T) = u^da and the brackets () around the 
indices to denote symmetrization, z.e., a^"b^^ = {a°'h^ + a^h°')/2. 

These constraints are simply the equations of fluid dynamics expanded to second order 
in the derivative expansion. The first few terms in the expansion ( |2.8| ) are familiar. The 
derivative free terms describe a perfect fluid with pressure [i.e., negative free energy den- 
sity) tt'^T'^, and so (via thermodynamics) entropy density s = An'^T^. The viscosity 77 of 



^"^Note that the variation in the radial direction, r, is never slow. Although we work order by order in 
the field theory derivatives, we will always solve all differential equations in the r direction exactly. 

^^Throughout this paper T'^" = IGnG^t'^'^ where G5 is the five dimensional Newton and t^'^ is the 
conventionally defined stress tensor, i.e., the charge conjugate to translations of the coordinate v. 
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this fluid may be read off from the coefficient of a^^'^ and is given by tt^T^. Notice that 
Tj/ s = l/(47r), in agreement with the famous result of Policastro, Son and Starinets IQ. 

Our computation of the two derivative terms in ( |2.8| ) is new; the coefficients of these 
terms are presumably related to the various 'relaxation times' discussed in the literature 
(see for instance [0). As promised earlier, the fact that we are dealing with a particular 
conformal fluid, one that is dual to gravitational dynamics in asymptotically AdS space- 
times, leads to the coefficients being determined as fixed numbers. It would be interesting 
to check whether the stress tensor determined above fits into the framework of the so called 



Israel-Stewart formalism ||39[ (see [^, ^ for reviews). R. Loganayagam is currently 
investigating this issue. 



In § |6.1| we have checked that the minimal covariantization of the stress tensor (|2.8| ) 
transforms as T^'^ — > e~^'^T^^ under the Weyl transformation 7]^^, e^'^^T^^j,, T — > e~'^T, 

— e~'^u'^, for an arbitrary function <^(x'^).[^ This transformation (together with the 
manifest tracelessness of T^'^) ensures Weyl invariance of the fluid dynamical equation 
( p.7|) . Note that we have computed the fluid dynamical stress tensor only in flat space. 
The generalization of our expression above to an arbitrary curved space could well in- 
clude contributions proportional to the spacetime curvature tensor. The fact that ( |2.8| ) 
is Weyl invariant by itself is a bit of a (pleasant) surprise. It implies that that the sum 
of all curvature dependent contributions to the stress tensor must be independently Weyl 
invariant. 



3 The perturbative expansion 

As we have described in § |], our goal is to set up a perturbative procedure to solve Ein- 
stein's equations in asymptotically AdS spacetimes order by order in a boundary derivative 
expansion. In this section we will explain the structure of this perturbative expansion, and 
outline our implementation of this expansion to second order, leaving the details of com- 
putation to future sections. 

3.1 The basic set up 

In order to mathematically implement our perturbation theory, it is useful to regard b and 
Pi described in § |] as functions of the rescaled field theory coordinates e where e is a 
formal parameter that will eventually be set to unity. Notice that every derivative of /3j 
or b produces a power of e, consequently powers of e count the number of derivatives. We 
now describe a procedure to solve Einstein's equations in a power series in e. Consider the 

R. Loganayagam gl) infor ms us that he has succeeded in rewriting our stress tensor in a number of 
different compact forms, one of which makes its covariance under Weyl transformations manifest. 
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metricp^ 

g = g^'\f3.,b)+eg^'\(3,,b)+e'g^'\(3.,b) + O {e^) , (3.1) 

where g^^^ is the metric (|2.6|) and g^^\g^'^^ etc are correction metrics that are yet to be 
determined. As we will explain below, perturbative solutions to the gravitational equations 
exist only when the velocity and temperature fields obey certain equations of motion. These 
equations are corrected order by order in the e expansion; this forces us to correct the 
velocity and temperature fields themselves, order by order in this expansion. Consequently 
we set 

= /5f ) + e + O {e^) , h = 6(0) + e fe^^) + O {e') , (3.2) 

where /Sf'"'' and 6'-") are all functions of ex'^. 

In order to proceed with the calculation, it will be useful to fix a gauge. We work with 
the 'background field' gauge 

grr = 0, 9r^^u^, Tr ((^7 ^ ) " V"^ = V n > 0. (3.3) 

Notice that the gauge condition at the point x'^ is given only once we know n^(w,x*). In 
other words, the choice above amounts to choosing different gauges for different solutions, 
and is conceptually similar to the background field gauge routinely used in effective action 
computations for non abelian gauge theories. 

3.2 General structure of perturbation theory 

Let us imagine that we have solved the perturbation theory to the (n — 1)**^ order, i.e., 
we have determined g^"^^ for m < n — 1, and have determined the functions Z?}"*^ and fc'-™-' 
for m < n — 2. Plugging the expansion ( p.l|) into Einstein's equations, and extracting the 
coefficient of e", we obtain an equation of the form 

g^^\xn = Sn. (3.4) 
Here if is a linear differential operator of second order in the variable r alone. As 

gin) is 

already of order e"', and since every boundary derivative appears with an additional power 
of e, H is an ultralocal operator in the field theory directions. It is important to note that 
i7 is a differential operator only in the variable r and does not depend on the variables x'^. 
Moreover, the precise form of this operator at the point x'^ depends only on the values of 
pf'^ and fe^O'' at but not on the derivatives of these functions at that point. Furthermore, 
the operator H is independent of n; we have the same homogeneous operator at every 

^•^For convenience of notation we are dropping the spacetime indices in 5^"^. We also suppress the 
dependence of b and /3,; on x^. 



H 



,(0) (^(0)^^(0)) 
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order in perturbation theory. 

The source term s„ however is different at different orders in perturbation theory. It 
is a local expression of n^^ order in boundary derivatives of (3f^ and 6^°^ as well as of 
{n — ky^ order in (3l''\ h^^^ for all /c < n — 1. Note that [3^^^ and 6'-"^ do not enter the ri"^ 
order equations as constant (derivative free) shifts of velocities and temperatures solve the 
Einstein's equations. 

The expressions (|3.4|) form a set of 5 x 6/2 = 15 equations. It turns out that four of 
these equations do not involve the unknown function g^"''^ at all; they simply constrain the 
velocity functions h and There is one redundancy among the remaining 11 equations 
which leaves 10 independent 'dynamical' equations. These may be used to solve for the 10 
unknown functions in our gauge fixed metric correction g^'^\ as we describe in more detail 
below. 

3.2.1 Constraint equations 

By abuse of nomenclature, we will refer to those of the Einstein's equations that are of 
first order in r derivatives as constraint equations. Constraint equations are obtained by 
dotting the tensor Emn with the vector dual to the one-form dr. Four of the five constraint 
equations (z.e., those whose free index is a /i index) have an especially simple boundary 
interpretation; they are simply the equations of boundary energy momentum conservation. 
In the context of our perturbative analysis, these equations simply reduce to 

^.^(ri) = (3.5) 

where T^'^_^-^ is the boundary stress tensor dual the solution expanded up to 0{e'^~^). 
Recall that each of g^^\g^^\.. are local functions of 6, /Sj. It follows that the stress tensor 
'^(n-i) ^Is*^ ^ local function (with at most n — 1 derivatives) of these temperature and 
velocity fields. Of course the stress tensor T^_^~^ also respects 4 dimensional conformal 
invariance. Consequently it is a 'fluid dynamical' stress tensor with n — 1 derivatives, the 
term simply being used for the most general stress tensor (with n—1 derivatives), written 
as a function of and T, that respects all boundary symmetries. 

Consequently, in order to solve the constraint equations at n^^ order one must solve 
the equations of fluid dynamics to (n — l)*'^ order. As we have already been handed a 
solution to fluid dynamics at order n — 2, all we need to do is to correct this solution to 
one higher order. Though the question of how one goes about improving this solution is 
not the topic of our paper (we wish only to establish a map between the solutions of fluid 
mechanics and gravity, not to investigate how to find the set of all such solutions) a few 
words in this connection may be in order. The only quantity in ( p.5|) that is not already 
known from the results of perturbation theory at lower orders are (3^"''^^ and b^'^~^\ The 
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four equations ( p.5| ) are linear differential equations in these unknowns that presumably 
always have a solution. There is a non-uniqueness in these solutions given by the zero 
modes obtained by linearizing the equations of stress energy conservation at zeroth order. 
These zero modes may always be absorbed into a redefinition of I3f'\ ft*-"-*, and so do not 
correspond to a physical non-uniqueness {i.e., this ambiguity goes away once you specify 
more clearly what your zeroth order solution really is). 

Our discussion so far may be summarized as follows: the first step in solving Einstein's 
equations at -n}^ order is to solve the constraint equations - this amounts to solving the 
equations of fluid dynamics at {n — 1)**^ order ( p.5|) . As we explain below, while it is of 
course difficult in general to solve these differential equations throughout R^'^, it is easy 
to solve them locally in a derivative expansion about any point; this is in fact sufficient to 
implement our ultralocal perturbative procedure. 



3.2.2 Dynamical equations 

The remaining constraint E„ and the 'dynamical' Einstein's equations E^^ may be used 
to solve for the unknown function g^'^\ Roughly speaking, it turns out to be possible to 
make a judicious choice of variables such that the operator H is converted into a decoupled 
system of first order differential operators. It is then simple to solve the equation (O) 



for an arbitrary source s„ by direct integration. This procedure actually yields a whole 
linear space of solutions. The undetermined constants of integration in this procedure are 
arbitrary functions of and multiply zero modes of the operator ( |3.4| ). As we will see 
below, for an arbitrary non-singular and appropriately normalizable source s„ (of the sort 
that one expects to be generated in perturbation theory]^), it is always possible to choose 
these constants to ensure that 

gin) is 

appropriately normalizable at r = oo and non-singular 
at all nonzero r. These requirements do not yet completely specify the solution for g^'^\ as 
H possesses a set of zero modes that satisfy both these requirements. A basis for the linear 
space of zero modes, denoted Qb and Qi, is obtained by differentiating the 4 parameter class 
of solutions ( p.2|) with respect to the parameters h and j3i. In other words these zero modes 
correspond exactly to infinitesimal shifts of jdf^ and 6^°-* and so may be absorbed into a 
redefinition of these quantities. They reflect only an ambiguity of convention, and may be 
fixed by a 'renormalization' prescription, as we will do below. 



I'^Provided the solution at order n — 1 is non-singular at all nonzero r, it is guaranteed to produce a 
non-singular source at all nonzero r. Consequently, the non-singularlity of s„ follows inductively. We think 
is possible to make a similar inductive argument for the large r behaviour of the source, but have not yet 
formulated this argument precisely enough to call it a proof. 
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3.2.3 Summary of the perturbation analysis 

In summary, it is always possible to find a physically unique solution for the metric g^'^\ 
which, in turn, yields the form of the n^^ order fiuid dynamical stress tensor (using the 
usual AdS/CFT dictionary). This process, being iterative, can be used to recover the fiuid 
dynamics stress tensor to any desired order in the derivative expansion. 

In § and § we will provide a few more details of our perturbative procedure, 
in the context of implementing this procedure to first and second order in the derivative 
expansion. 

3.3 Outline of the first order computation 

We now present the strategy to implement the general procedure discussed above to first 
order in the derivative expansion. 



3.3.1 Solving the constraint equations 

The Einstein constraint equations at first order require that the zero order velocity and 
temperature fields obey the equations of perfect fiuid dynamics 



S^Tj; = , (3.6) 

where up to an overall constant 

^(o; = p)jl(^'^^ + 4<o)«(o))- (3.7) 

While it is difiicult to find the general solution to these equations at all x^, in order to 
carry out our ultralocal perturbative procedure at a given point y^, we only need to solve 
these constraints to first order in a Taylor expansion of the fields b and (3i about the point 
y^. This is, of course, easily achieved. The four equations ( p.6|) may be used to solve for 
the 4 derivatives of the temperature field at y'^ in terms of first derivatives of the velocity 
fields at the same point. This determines the Taylor expansion of b to first order about y'^ 
in terms of the expansion, to first order, of the field /5j about the same point. We will only 
require the first order terms in the Taylor expansion of velocity and temperature fields in 
order to compute g^^^y^)- 



3.3.2 Solving the dynamical equations 

As described in the previous section, we expand Einstein's equations to first order and 
find the equations (|3^ ). Using the 'solution' of § |3.3.1| , all source terms may be regarded 
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as functions of first derivatives of velocity fields only. The equations ( |3.4| ) are then easily 
integrated subject to boundary conditions and we find ( |3.4| ) is given by 



fif^^^ = Qp^ + fb{Xi, V) Qh + fi{Xj, V) Qi, 



{3.i 



where Qp^ is a particular solution to ( |3.4D , and and fi are a basis for the zero modes of 
H that were described in the § |3.2|. Plugging in this solution, the full metric g^^^ + g^^\ 



when expanded to order first order me, is (|3.1|) 

9 = 9^'^ + e (9? + {h + h^''')9i> + (/. + (^^'^)9^ , , 



(3.9) 



where the four functions of x'^, fb + fc*-^'', fi + (3^^^ are all completely unconstrained by the 
equations at order e. 



3.3.3 The 'Landau' Frame 



Our solution ( |3.8|) for the first order metric has a four function non-uniqueness in it. As 
fb and fi may be absorbed into b^^^ and (3^^^ this non-uniqueness simply represents an 
ambiguity of convention, and may be fixed by a 'renormalization' choice. We describe our 
choice below. 

, it is straightforward to use the AdS/CFT correspondence to recover the 
To first order in e the boundary stress tensor dual to the metric (|3!^ 



Given g^^'' 
stress tensor, 
evaluates to 



rpfll/ 



where 



l(r,- + 4«'^«^)-|T-, 



6 = 6(°)+£(6(^) + /b) 

A = ^(/Sf ^ + f^) 



(3.10) 



(3.11) 



where Tj^^^, defined by ( 3.10| ), is an expression linear in derivatives of the velocity fields 
and temperature fields. Notice that our definition of T^/^^, via (|3.10|) , depends explicitly on 
the value of the coefficients fi, ft of the homogeneous modes of the differential equation 
( p.4|) . These coefficients depend on the specific choice of the particular solution gp\ which is 
of course ambiguous up to addition of homogenous solutions. Any given solution ( |3.8|) may 
be broken up in many different ways into particular and homogeneous solutions, resulting 
in an ambiguity of shifts of the coefficients of fb, fi and thereby an ambiguity in Tj/^^. It 
is always possible to use the freedom provided by this ambiguity to set U(^o)^Tl^^ = 0. 

This choice completely fixes the particular solution gp\ We adopt this convention for the 
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particular solution and then simply simply set g''^^ = Qp ^-c, choose fb = fi = 0. Tj^^ is 
now unambiguously defined and may be evaluated by explicit computation; it turns out 
that 

i(^) - a . 

The discussion of the previous paragraph has a natural generalization to perturbation 
theory at any order. As the operator H is the same at every order in perturbation theory, 
the ambiguity for the solution of (^'•"^ in perturbation theory is always of the form described 
in ( p.9|) . We will always fix the ambiguity in this solution by choosing u^T^^^-^ = 0. The 
convention dependence of this procedure has a well known counterpart in fiuid dynamics; 
it is simply the ambiguity of the stress tensor under field redefinitions of the temperature 
and u^. Indeed this field redefinition ambiguity is standardly fixed by precisely the 'gauge' 
choice u^Tl^^ = 0. This is the so called 'Landau frame' widely used in studies of fiuid 
dynamics.0 

We present the details of the first order computation in § |^ below. 



3.4 Outline of the second order computation 

Assuming that we have implemented the first order calculation described in § |3.3| , it is 
then possible to find a solution to Einstein's equations at the next order. In this case care 
should be taken in implementing the constraints as we discuss below. 



3.4.1 The constraints at second order 



The general discussion of § ^]2| allows us to obtain the second order solution to Einstein's 



equations once we have solved the first order system as outlined in § |3.3| . However, we 
need to confront an important issue before proceeding, owing to the way we have set up 
the perturbation expansion. Of course perturbation theory at second order is well defined 
only once the first order equations have been solved. While in principle we should solve 
these equations everywhere in R^'^, in the previous subsection we did not quite achieve 
that; we were content to solve the constraint equation (p.6|) only to first order in the Taylor 
expansion about our special point y^. While that was good enough to obtain g^^\ in order 
to carry out the second order calculation we first need to do better; we must ensure that 
the first order constraint is obeyed to second order in the Taylor expansion of the fields 
and ^ about y^. That IS, we require 

9a9^TJ) (2/") = 0- (3-12) 



^^Conventionally, one writes in fluid mechanics the stress tensor as the perfect fluid part and a dissipative 
part i.e., T'"" = T^^j^^^^ + Tj^^^^.^^^^^^. The Landau gauge condition we choose at every order simply 
amounts to u,, T^"' ^- =0. 

A* dissipative 
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Essentially, we require that T^-q^ satisfy the conservation equation ( |3.6| ) to order before 
we attempt to find the second order stress tensor. In general, we would have need ( |3.6| ) 
to be satisfied globally before proceeding; however, the ultralocality manifest in our set-up 
implies that it suffices that the conservation holds only to the order we are working. If we 
were interested in say the n*^ order stress tensor Tl^^-^ we would need to ensure that the 
stress tensor up to order n — 1 satisfies the conservation equation to order O 

The equations (|3.12|) may be thought of as a set of 16 linear constraints on the coeffi- 
cients of the (40-1-78) two derivative terms involving 6^°^ and (3f'\ We use these equations 
to solve for 16 coefficients, and treat the remaining coefficients as independent. This pro- 
cess is the conceptual analogue of our zeroth order 'solution' of fluid dynamics at the point 

(described in the previous subsection), obtained by solving for the flrst derivatives of 
temperature in terms of the flrst derivatives of velocities. Indeed it is an extension of that 
procedure to the next order in derivatives. See § |^ for the details of the implementation 
of this procedure. In summary, before we even start trying to solve for g^'^\ we need to 
plug a solution of ( p.l2|) into g^^^ + g^^^ expanded in a Taylor series expansion about y^. 
Otherwise we would be expanding the second order equations about a background that 
does not solve the flrst order fluid dynamics. 



3.4.2 Nature of source terms 

As we have explained above, the Einstein's equations, to second order, take the schematic 
form described in ( ^.4| ) 

9^'^=Sa + Sb (3.13) 

We have broken up the source term above into two pieces, Sa and Sb, for conceptual 
convenience. Sa is a local functional of f3^^^ and fe^"-* of up to second order in fleld theory 
derivatives. Terms contributing to Sa have their origin both in two fleld theory derivatives 
acting on the metric g^^^ and exactly one fleld theory derivative acting on g^^^ (recall that 
g^^^ itself is a local function of jSf'^ and 6^°^ of flrst order in derivatives). The source term 
Sf, is new: it arises from flrst order derivatives of the velocity and temperature corrections 
Z?}^'* and b^^\ This has no analogue in the flrst order computation. 

As we have explained above, b^^^ are absolutely any functions that obey the equa- 
tions ( |3.6| ). In particular, if it turns of that the functions pf'^ + e jS^^^ and b^^'' + e fe*'^-' obey 
that equation (to flrst order in e) then f^l^'' and 6^^^ may each simply be set to zero by an 
appropriate redeflnition of /?|°^ and bf'\ This results in a 'gauge' ambiguity of the func- 
tions In our ultralocal perturbative procedure, we choose to flx this ambiguity 
by setting (3^^^ to zero (at our distinguished point y^) while leaving fe^^-' arbitrary.|3 

The functions 6^^\/3,p'' have sixteen independent first derivatives. aU but four of which may be fixed 
by the gauge freedom. We choose use this freedom to set ah velocity derivatives to zero. 



H 



,{0)(^(0)^^(0)) 
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3.4.3 Solution of the constraint equations 

With the source terms in place, the procedure to solve for g^'^^ proceeds in direct imitation 
of the first order calculation. The constraint equations reduce to the expansion to order e 
of the equation of conservation of the stress tensor 

T^-" = ^ (4 u^u" + 7]"") -2e^^a^"' (3. 14) 

with = (3^^\ h = +eb'-^\ These four equations may be used to solve for the four 
derivatives d^b^^^ at x^. Consequently the constraint equations plus our choice of gauge, 
uniquely determined the first order correction of the temperature field b^^^ and velocity 
field ^ as a function of the zeroth order solution. 

Note that the gauge (5^^\y^) = may be consistently chosen at any one point y^, 
but not at all x^. Nonetheless the results for g^'^^ that we obtain using this gauge will, 
when appropriately covariantized be simultaneously applicable to every spacetime point 
x'^. The reason for this is that all source terms depend on 6^^^ and (3^^^ only through the 
expansion to order e of d^T^'^ = with T^'^ given by ( gmi) . Note that this source term is 
'gauge invariant' (recall that 'gauge' transformations are simply shifts of b^^^ and (3^^^ by 
zero modes of this equation). It follows that g^"^^ determined via this procedure does not 
depend on our choice of gauge, which was made purely for convenience. 

3.4.4 Solving for g^'^^ and the second order stress tensor 

Now plugging this solution for b^^^ into the source terms it is straightforward to integrate 
( p.l3| ) to obtain g^'^\ We fix the ambiguity in the choice of homogeneous mode in this 
solution as before, by requiring T^g)' ""(0)1^ = 0. This condition yields a unique solution for 
g^"^^ as well as for the second order correction to the fluid dynamical stress tensor T^^, 
giving rise to the result (|2.9|). We present the details of the second order computation in 
§|. 

In the rest of this paper we will present our implementation of our perturbative proce- 
dure described above, to first and second order in the derivative expansion. 

4 The metric and stress tensor at first order 

In this section we will determine the solution, to first order in the derivative expansion. As 
we have described in § ^, the equations that determine g^^^ at x^ are ultralocal; consequently 
we are able to solve the problem point by point. It is always possible to choose coordinates 
to set = (1,0,0,0) and 6^°) = 1 at any given point x^. Making that choice, the metric 
( p.6|) expanded to first order in derivatives in the neighbourhood of (chosen to be the 
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origin of R^'^ for notational simplicity) is given by 



ds f^Q^^ = 2 dv dr — f [r] dv +r dxidx^ 

T-M (9 a(o) (4-1) 
- 2 d^(3f^ dx' dr-2 x^'d^lSf^ r^ (1 - /(r)) dx' dv - 4 — -| dv^ . 

In order to implement the perturbation programme described in the previous section, we 
need to find the first order metric g^^^ which, when added to ( [4. 11 ), gives a solution to 
Einstein's equations to first order in derivatives. 

The metric ( [4.1|) together with g^^'' has a background piece (the first line in (f4.1|) ) 



which is simply the metric of a uniform black brane. In addition it has small first derivative 
corrections, some of which are known (the second line of ( |4.1D ), and the remainder of which 
((yf*^^^) we have to determine. Now note that the background black brane metric preserves 
a spatial 5*0(3) rotational symmetry. This symmetry allows us to solve separately for the 
50(3) scalars, the 5*0(3) vector and 50(3) symmetric traceless two tensor (5) components 
of g^^^ and lies at the heart of the separability of the matrix valued linear operator H into 
a set of ordinary linear operators. 

In the following we will discuss each of these sectors separately and determine g^^\ 



Subsequently, in § 4.4 we present the full solution to order e and proceed to calculate the 



stress tensor in 5 H75 



4.1 Scalars of 50(3) 

The scalar components of g^^^ are parameterized by the functions hi(r) and ki{r) according 



^i^^(^) = Sr^ hi{r) 

9iv{r) = ^ (4.2) 



2 



Here g^^ and related to each other by our gauge choice Tr((^(°))-^^(^)) = 0. 

The scalar Einstein's equations {i.e., those equations that transform as a scalar of 
50(3)) may be divided up into constraints and dynamical equations. The constraint 
equations are obtained by contracting Einstein's equations (the first line of ( ^.1| )) with the 



^^In the spatial R'^ C R'^'^ we will often for ease of notation, avoid the use of covariant and contravariant 
indices and adopt a summation convention for repeated indices i.e., g^P = 9iP ■ 
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vector dual to the one form dr. The first scalar constraint is 



r2 /(r) E,r + E,, = 0, (4.3) 

which evaluates to 

= ^ . (4.4) 

Below, we will interpret ( [4.4|) as the expansion of the fluid dynamical stress energy conser- 
vation, expanded to first order. The second constraint equation, 

r^f{r)Err + E,r = 0, (4.5) 

leads to 

12r3/ii(r) + (3r^ - 1) h[{r) - k[{r) = -Qr^ . (4.6) 

3 

To this set of constraints we need add only one dynamical scalar equation,Q the simplest 
of which turns out to be 

5h[{r) + rh';{r) = . (4.7) 

The LHS of ([4.7|) and ([4.6|) are the restriction of the operator H of ( p. 41) to the scalar 
sector. The RHS of the same equations are the scalar parts of the source terms Si. Notice 
that if is a first order operator in the variables h[{r) and ki{r). Consequently the equation 
( [4.7|) may be integrated for an arbitrary source term. The resulting solution is regular at all 
nonzero r provided that the source shares this property, and the growth hi{r) at infinity 
is slower than a constant - the behaviour of a non normalizable operator deformation 
- provided the source in ([4.7|) grows slower than 1/r at large r. Once hi{r) has been 
obtained ki{r) may be determined from ([4.6| ) by integration, for an arbitrary source term. 
Once again, the solution will be regular and grows no faster than at large r, provided 
the source in that equation is regular and normalizable. The two source terms of this 
subsection satisfy these regularity and growth requirements, and it seems clear that this 
result will extend to arbitrary order in perturbation theory (see the next section). 

The general solution to the system ( |4.(j| ) and ( [4.7| ), obtained by the integration described 
above, is 



+ 3r S- — 

3 



hi{r) = s + —, ki{r) = -f^ + 3r\s - — + u , (4.8) 



where s,t and u are arbitrary constants (in the variable r). In the solution above, the 
parameter s multiplies a non normalizable mode (which represents a deformation of the field 



""^^We have explicitly checked that the equations listed here imply that the second dynamical equation 
is automatically satisfied. 
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theory metric) and so is forced to zero by our boundary conditions. A linear combination of 
the pieces multiphed by t and u is generated by the action of the coordinate transformation 
r' = r (1 + a/r^) and so is pure gauge, and may be set to zero without loss of generality. 
The remaining coefficient u corresponds to an infinitesimal temperature variation, and is 
forced to be zero by our renormalization condition on the stress tensor m^^^ T^,^ = (see 
the subsection on the stress tensor below). In summary, each of s,t,u may be set to zero 
and the scalar part of the metric g^^\ denoted g^g \ is 



Two comments about this solution are in order. First note that ki{r) is manifestly 
regular at the unperturbed 'horizon' r = 1, as we require. Second, it grows at large 
r like r^. This is intermediate between the r° growth of finite energy fluctuations and 
the growth of a field theory metric deformation. As g^^^ + g^^^ obeys the Einstein's 
equations to leading order in derivatives, the usual Fefferman- Graham expansion assures 
us that the sum of first order fluctuations in g^^^ + g^^^ must (in the appropriate coordinate 
system) die off like compared to terms that appear in the zeroth order metric (this 
would correspond to ki{r) constant at infinity). Consequently the unusually slow fall off 
at infinity of our metric g^^^ must be compensated for by an equal but opposite effect from 
a first order fluctuation piece in the second line of (|4.1| ). This indeed turns out to be the 
case. While an explicit computation of the boundary stress tensor dual to ( |4.1|) yields a 
result that diverges like r^, this divergence is precisely cancelled when we add g^^^ above 
to the metric, and the correct value of the stress dual to g^'^'' + g^^^ is in fact zero in the 
scalar sector, in agreement with our renormalization condition M(o)^t T'^'^ = 0. 

4.2 Vectors of SO (3) 

In the vector channel the relevant Einstein's equations are the constraint /(r) Eri + E^i = 
and a dynamical equation which can be chosen to be any linear combination of the 
Einstein's equations E^ = and Eyi = 0. The constraint evaluates to 



which we will later interpret as a consequence of the conservation of boundary momentum. 
In order to explore the content of the dynamical equation (we choose Eri = 0), it is 
convenient to parameterize the vector part of the fluctuation metric by the functions 

as 




(4.9) 



(4.10) 




(4.11) 
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The dynamical equation for ji(r) turns out to be 



d / 1 d 



dr \ dr^ 



(1)/ 



—0 dv/3l 



(0) 



(4.12) 



The LHS of ( [4.12 ) is the restriction of the operator H of ( |3.4|) to the vector sector, and 
the RHS of this equation is the projection of si to the vector sector. H is of first order in 
the variable j^^^' (r) and so may be integrated for an arbitrary source term. The resulting 
solution is regular and normalizable provided the source is regular and decays at infinity 
faster than 1/r. This condition is obeyed in ( [4.12| ); it seems rather clear that it will continue 
to be obeyed at arbitrary order in perturbation theory (see the next section). 



Returning to ([4.12|) , the general solution of this equation is 



(4.13) 



for arbitrary constants Oj, q. The coefficient Oj multiplies a non-normalizable metric defor- 
mation, and so is forced to zero by our choice of boundary conditions. The other integration 
constant q multiplies an infinitesimal shift in the velocity of the brane. It turns out (see 
below) that a nonzero value for Cj leads to a nonzero value for Toj which violates our 
'renormalization' condition, consequently q must be set to zero. In summary, 

(^g'^^^ ^ dx^dx^ = 2 r ^ dv dx\ (4.14) 

As in the scalar sector above, this solution grows by a factor of faster at the boundary 
than the shear zero mode. This slow fall off leads to a divergent contribution to the stress 
tensor which precisely cancels an equal and opposite divergence from terms in the expansion 
of g^^^ to first order in derivatives. As we will see below, the full contribution of (7^°^ + g^^^ 
to the vector part of the boundary stress tensor is just zero, again in agreement with our 
renormalization conditions. 



4.3 The symmetric tensors of SO (3) 

We now turn to g!^\ the part of g^^^ that transforms in the 5, the symmetric traceless two 
tensor representation, of SO (3). Let us parameterize our metric fiuctuation by 

(^g?^ ^ rfx" dxf^ = ajj^ (r) dx' dx\ (4.15) 
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where aij is traceless and symmetric. The Einstein's equation Eij = yield 



± (r^ fir) ±a^^) 



R 2 (0) 

-or a- 



(4.16) 



where we have defined a symmetric traceless matrix 



cr, 



(0) 



(4.17) 



The LHS of ( 4.16| ) is the restriction of the operator H of ( |3.4|) to the tensor sector, and 
the RHS of this equation is the tensor part of the source term s\. Note that if is a first 
order operator in the variable a 



and so may be integrated for an arbitrary source 
term. The solution to this equation with arbitrary source term s(r) is given by (dropping 
the tensor indices): 

Note that the lower limit of the inner integral in ([4.18|) has been chosen to be unity. 
Provided that s[x) is regular at x = 1 (this is true of ( [4.16| ) and will be true at every order 
in perturbation theory), s{x) has a zero at x = 1. It follows that the outer integrand 
in ( [4.18[ ) is regular at nonzero x (and in particular at x = 1) despite the explicit zero in 
the factor f{x) in the denominator. The solution for is also normalizable provided the 
source is regular and grows at infinity slower than r^. This condition is obeyed in ( [4. 161 ) 
and is expected to continue to be obeyed at arbitrary order in perturbation theory (see the 
next section). 



Applying ( [4. 18] ) to the source term in ( |4.16| ) we find that the solution for a^-j^ is given 



by 



(9^ 



with 



F(r) 



dx 



X^ + X + 1 



x{x + 1) (a;2 + 1) 



dx°'dx' 



1 
4 



2r'^F{r)af dx'dx^. 



In 



(l + r)2(l + r2) 



— 2 arctan(r) + vr 



(4.19) 



(4.20) 



At large r it evaluates to 



{g'>P)ai3 dx'^dx'^ 



(tJ°^ dx'dx^. 

''J 



(4.21) 



As in the previous subsections, the first term in ( [4.21|) yields a contribution to the stress 
tensor that diverges like r^, but precisely cancels the corresponding divergence from first 
derivative terms in the expansion of g^^\ However the second term in this expansion yields 
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an important finite contribution to the stress tensor, as we will see below. 



Summary of the first order calculation: In summary, our final answer for g^^^ + g^^\ 
expanded to first order in boundary derivatives about = 0, is given explicitly as 

ds^ = 2dv dr — r'^f{r) dv"^ + dxi dx^ 

-2x^ d^.pf'^ dr dx' -2x^df, /9j°V2(l - /(r)) dv dx' - 4 — — dv^ (4.22) 

+ 2 F(r) ajf dx' dx^ + ^ r dipf'^ dv^ + 2r d^fif^ dv dx\ 

This metric solves Einstein's equations to first order in the neighbourhood of = 
provided the functions h^^^ and j3f''^ satisfy 

Q 5(0) _ <=^^P^ 

3 (4.23) 
4.4 Global solution to first order in derivatives 

In the previous subsection we have computed the metric g^^^ about x^ assuming that 
6(0) = 1 and = at the origin. Since it is possible to choose coordinates to set an 
arbitrary velocity to zero and an arbitrary fe'-O) to unity at any given point (and since our 
perturbation procedure is ultralocal), the results of the previous subsection contain enough 
information to write down the metric g'^^^ about any point. A simple way to do this is to 
construct a covariant metricf^, as a function of and 6, which reduces to ([4.22|) when 
6'-°) = 1 and (3^^ = 0. It is easy to check that 

ds'^ = —2 dx'^dr — f{br) u^u^ dx'^dx" + P^^ dx'^dx" 

2 (4.24) 

+ 2r'^hF{hr) a^^ dx^dx" + - r u^u^ dxu'^ dx'^dx" - r u^dx (u^u^) dx^dx", 

3 

does the job, up to terms of second or higher order in derivatives. Here we have written 
the metric in terms of a^^^, defined in ( |2.9| ) and the function F{r) introduced in ( |4.20| ). 
Furthermore, it is easy to check that the metric above is the unique choice respecting the 
symmetries (again up to terms of second or higher order in derivatives). It follows that 
( [4.24| ) is the metric g^'^^ + g^^'^- It is also easily verified that the covariant version of ( |4.23 ) 



^^By abuse of notation, we will refer to expressions transformation covariantly in the boundary metric 
(chosen here to be rj^^) as covariant. In particular, we are not interested in full bulk covariance as we will 
continue to restrict attention to a specific coordinatization of the fifth direction. 



22 



is ( |3.6| ). We will interpret this as an equation of stress energy conservation in the next 
subsection. 



4.5 Stress tensor to first order 

Given the solution to the first order equations, we can utilize the AdS/CFT dictionary to 
construct the boundary stress tensor using the prescription of |^ (see also ^^). For the 
metric ([4 .241) it is not difficult to compute the stress tensor; all we need to do is compute 
the extrinsic curvature tensor K^j^^y to the surface at fixed r. By convention, we choose the 
unit normal to this surface to be outward pointing, i.e. pointing towards the boundary, in 
the definition of K^^. Using then the definition 

T^' = -2 lim (K^ - 6^) , (4.25) 



+ 00 



on our solution ([4.24|) , we find the result is given simply as 



T'"' = (4 u'^u" + 7]^ - ^ a^r (4.26) 



where a^'^ was defined in ( p.9| ) and all field theory indices are raised and lowered with the 
boundary metric rj^j^^. As explained in the introduction, this stress tensor implies that the 
ratio of viscosity to entropy density of our fluid is l/(47r). Note that as mentioned previ- 
ously, the expression ( [4.26] ) is only correct up to first derivative terms in the temperature 



{T = 1/b) and velocities. 



5 The metric and stress tensor at second order 

In order to obtain the metric and stress tensor at second order in the derivative expansion, 
we follow the method outlined in § ^ and implemented in detail in § ^ to leading order. 
Concretely, we choose coordinates such that f^f^^ = and b^^^ = 1 at the point = 0. 
The metric g^^^ + g^^^ given in ( |4.24| ) may be expanded to second order in derivatives. 



This involves Taylor expanding g^^^ to second order and g^^^ to first order, the second order 
analogue of ( [4.1[ ). As we have explained in § ^.4| , at this stage we also make the substitution 
^(0) jjio) _|_ 5(1) ^ a^]2d treat b^^^ as an order e term, and so retain only those expressions 
that are of first derivative order in b^^^ (and contain no other derivatives). This process is 
straightforward and we will not record the (rather lengthy) resultant expression here. To 
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this expression we add the as yet undetermined metric fluctuation 



dx'^dx^ = -3 /i2(r) dv dr + /lafr) dxi dx' + dv^ + 2 dv dx' + r'^af^ dx' dxK 

(5.1) 



We plug this metric into Einstein's equations and obtain a set of linear second order 
differential equations that determine h2,k2,jl yC^lj ■ As in the previous section, 50(3) 

(2) 

symmetry ensures that the equations for the scalars /i2, k2, the vectors jl , and the tensor 
a'lj^ do not mix. Moreover, as we have explained in § the equations that determine 
these unknown functions are identical to their first order counterparts in the homogeneous 
terms, but differ from those equations in the sources. As a result, the only new calculation 
we have to perform in order to obtain the metric at second order is the computation of the 
source terms. Once these terms are available, the corresponding equations may easily be 
integrated, as in the previous section. 

Recall that the input metric into Einstein's equations includes terms that arise out 
of the Taylor expansion of g^'^^ + g^^^ that have explicit factors of the coordinates x^. 
Nonetheless, a very simple argument assures us that the source terms in the equations 
that determine g^"^^ must all be independent of x^. The argument runs as follows: We 
have exphcitly constructed g^^^ in the previous section so that Emn {g^^^ +9^^^) = Omn 
where Omn is a local expression constructed out of second order or higher x^ derivatives 
of velocity and temperature fields. It follows that x^ dependence of sources, which may be 
obtained by Taylor expanding Omn about x'^ = 0, occurs only at the three derivative level 
or higher. It follows that source terms at the two derivative level have no x'^ dependence. 
Clearly, this argument has a direct analogue at arbitrary order in perturbation theory. 

A crucial input into the argument of the last paragraph was the fact that g^^^ + g^^^ 
satisfies Einstein's equations in a neighbourhood of = (and not just at that point). 
As we have seen in the previous section, the fact that the energy conservation equation is 
obeyed ai x^ = allows us to express all first derivatives of temperature in terms of first 
derivatives of velocities (see ([4.10|) and ([4.4|) ). In addition, (3^^^ and b^^^ must be chosen 



so that (|3.12| ) is satisfied. The sixteen equations (|3.12|) can be grouped into sets that 



transform under 5*0(3) as two scalars, three vectors and one tensor [i.e., 5). We will now 
explain how these constraints may be used to solve for 16 of the independent expressions 
of second order in derivatives of velocity and temperature fields. 

In order to do this, let us first list all two derivative 'source' terms that can be built 
out of second derivatives of b^^^ or or out of squares of first derivatives of These 
expressions may be separated according to their transformation properties under 50(3) as 
scalars, vectors and tensors and higher order terms. The higher order pieces will not be 
of interest to us. An exhaustive list of these expressions that transform in the 1, 3 or 5 is 
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given in Table l.f^ We define the vector ii as the curl of the velocity i.e., 



1 of so {3) 


3 of 50(3) 


5 of 50(3) 


si = ld'^b 


vii = I didyb 


tiij = \didjh - \ ^3^i3 


S2 = dydiPi 






S3 = 1 d% 


V3i = dJi 


^3*^ dyCTij 


&i = d^PidyPi 


\4i = Idjaji - d'^pi 


liij = dy(3i dy(3j - 1 ©1 


&2 = ii dypi 


V5i = d'^pi 


%2ij = i{i dy(3j) - 1 ©2 5ij 


©o = (diPi? 




13,-,- = 2 e,.,r,- d.S'^ d^\B^ + 1 ©2 


©4 = ii f 










~ 3 ®4 ^ij 






"ZQij = (Jik (y^j - 1 ©5 ^ij 









Table 1: An exhaustive list of two derivative terms in made up from the temperature 
and velocity fields. In order to present the results economically, we have dropped the 
superscript on the velocities /3j and the inverse temperature b, leaving it implicit that these 
expressions are only valid at second order in the derivative expansion. 



i, = e,jkd^P'' , (5.2) 

and the symmetric traceless tensor (Tij has been previously defined in (|4.17|) . 

As a simple check on the completeness of expressions in Table 1, notice that the number 
of degrees of freedom in those of the tabulated expressions that are formed from a product 
of two single derivatives is 5 (in the scalar sector), 5 x 3 (in the vector sector), and 7 

^"^Note that the tensors are symmetric in their indices. The symmetrization as usual is indicated by 
parentheses. 
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X 5 in the tensor sector, leading to a total of 55 real parameters. Together with degrees 
of freedom from the two 7s and one 9 that can also be formed from the product of two 
derivatives (but will play no role in our analysis) this gives 78 degrees of freedom. This 
is in agreement with the expected | x 12 x 13 =78 ways of getting a symmetric object 
from twelve parameters (the first derivatives of the velocity fields) . On the other hand, the 
genuinely two derivative terms in Table 1 have 3x1 + 5x3 + 3x5 = 33 degrees of freedom 
which together with a two derivative term that transforms in the 7 (which however plays 
no role in our analysis) is the expected number 40 = 10 x 4 of two derivative terms arising 
from temperature and velocity fields. 

Assuming that we have already employed the first order conservation equation ( p. 61) to 
eliminate the first derivatives of b, we have to deal with the constraint equation (|3.12|) at 
the second order. Using the list of second order quantities given in Table 1, it is possible to 
show that (|3.12|) take the form of the following linear relations between these two derivative 
terms: 

51 = o ^3 - ©1 + n ®3 + ^ ©4 - o ®5 

3 9 6 3 

52 = S3 - 6i + - 64 - 65 

10 1 1^ 1^ 2^ 

vii = Y + - V5i + - Q3i, - - 5J2i - - QJ3i 

10 1 2^ 1 5^ 

V2i = — + - V5i - - + - 2J2i - - QJ3i 
9 9 3 6 3 

V3i = -^^4i + ^5i 



(5.3) 



Given these relations we now proceed to analyze the potential source terms arising from 
the metric ( [4.24| ) at (^(e^). The analysis, as before, can be done sector by sector - the 



computations for the scalar, vector and tensor sectors are given in § |5.1| , § |5]^ and § p!3| 
respectively. 
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5.1 Solution in the scalar sector 



Given the general second order fluctuation ( |5.1|) , we parameterize scalar components of 
g^"^^ in terms of the functions /i2(r) and k2{r) according to 

9S\r) = Sr^ /i2(r) 

S?M = (5.4) 
9S>(.r) = -|A2(r) . 



As we have explained in the § iA, the constraint Einstein's equations in this sector are 
given by the r and v component of the one-form formed by contracting the Einstein tensor 
with the vector dual to the one-form dr. The v component of this constraint, i.e. the 
second order expansion of ( |4.3| ), evaluates to 

jiya.i.") = ^65. (5.5) 

This equation enables us to solve for the flrst v derivative of b^^^ in terms of two derivative 
terms made up of f]f'\ but imposes no further constraints on b^^\ jS^^^ (|5.5|) has a simple 
physical interpretation; it is simply the time component of the conservation equation for 
the stress tensor ( [4.26| ), expanded to second order in derivatives. Consequently ( |5.5| ) is the 



Navier Stokes equation! 

The r component of the constraint, i.e. ( [4 .51) , gives us one relation between the functions 
h2{r) and k2{r) and their derivatives. As in § |4.1|, to this constraint we must add one 
dynamical equation. We obtain the following equations 

5h',{r)+rh'i{r) = Sh{r) 

Kir) = Sk{r) (5.6) 



where 



^ 4r 2r 1 + 2 1 

S,{r) = -—ss + 2r&i-—Gs + ^-^&4+-Wk{r)e5 . 



(5.7) 
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The functions Wh{r) and Wk{r) are given by 



, , A (r^ + r + 1? -2(3r^ + 2r + 1) F(r 



W,(r 



3 r{r + lY{r^ + iy 

2 4 (r^ + r + 1) (3 r"^ - 1) F(r) - {2r^ + 2r^ + 2r3 - r 

3 r (r + 1) (r2 + l) 



As advertised, it is clear that the differential operator acting on the functions h2{r) and 
k2{r) is identical to the one encountered in the first order computation in § |4.1| . The 
equation (|5.6| ) can be explicitly integrated; to do so it is useful to record the leading large 
r behaviour of the source term Sh{r): 

(^64+^65) • (5.8) 



The first equation in ( |5.(j| ) can be integrated given this asymptotic value to obtain the 
leading behaviour of the function h2{r). One finds 

= -iTi + / \^ dy [S,{y) - - S^^ . (5.9) 

The integral expression above can be shown to be of O (r^^) and hence the asymptotic 
behaviour of h2{r) is controlled by Sh- Given h2{r), one can integrate up the second 
equation of ( |5.6|) for k2{r). The leading large r behaviour of the source term Sk{r) is given 
by 

Sk{r)^rS^ = r (-U3 + 2ei-l&3-l&4+l&5] , (5-10) 



9 " 6 



and hence we have 



^2 



Hr) = ^S^~ I dx {S,{x)-xSr) . (5.11) 



In this case the integral makes a subleading contribution starting at O {r As in § ^ 



we have chosen the coefficients of homogeneous solutions to this differential equation so as 
to ensure normalizability and vanishing scalar contribution to the stress tensor (according 
to our renormalization conditions). 
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5.2 Solution in the vector sector 



The analysis of the vector fluctuations at second order mimics the computation described 
in § The vector fluctuation in g^"^^ is chosen as described in ( ^.11) to be 

.-(2) 

^i? = • (5-12) 

Once again, the analysis is easily done by looking at the constraint equations which are 
obtained by contracting the tensor Emn with the vector dual to dr. The i^^ constraint 
equation evaluates to 

15 33 

18 dib^'^ = 5 V4i + 5 + 15 5Ji, - — 2J2i - y . (5.13) 

This equation allows us to solve for the spatial derivatives of b^^^ in terms of derivatives of 
/3f ^ and b^°\ (ISTTl is simply the expansion to second order in derivatives of the conserva- 



tion of momentum of the stress tensor ( [4.26| ). 

To complete the solution in the vector channel, we need to solve for j^'^\r), which can 
be shown to satisfy a dynamical equation 

d f I d .(2) 



,,V.3*^rWJ=B,M. (5.14) 

Note that the LHS of this expression has the vector part of the operator H acting on 
Here Bj(r) is the source term which is built out of the second derivative terms transforming 
in the 3 of 5*0(3) given in Table 1. 

B(r) - + (5 15) 



with 



= 4 (10 V4 + V5 + 3 5Ji - 3 532 - 6 5J3) 
B*^'^ = 9 (20V4-5 2J2-6 2J3), 



(5.16) 



and we have introduced the polynomial: 

p(r) =2r^ + 2r2 + 2r-3 . (5.17) 
Clearly p(r) determines the large r behaviour of the vector perturbation; asymptotically 
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B(r) — s> g^B°°. Hence, integrating ( |5.14| ) we find that j^'^\r) is given as 



2 rocj POO / -1 



TDOO 

3^' 



(5.18) 



where once again we have chosen the coefficients of homogeneous modes in order to main- 
tain normahzability and our renormalization condition. As with the first order computation 
described in § [4.2| , the solution ( [5. 181 ) makes no contribution to the stress tensor of the field 
theory. 



5.3 Solution in the tensor sector 



Finally, we turn to the tensor modes at second order where we shall recover the explicit 
form of the second order contributions to the stress tensor. Our task is now to determine 
the functions alf{r) in (|5.1| ). As in § [4.3| , in the symmetric traceless sector of 5*0(3) one 
has only the dynamical equation given by 



1 d 
2r dr 



r 



4 ) y') 



Aii(r) 



(5.19) 



where 



A.,(r) = a,(r) (Ti., + -T4., + t3., 



with the coefficient functions 



ai(r 



a5(r) 



3p{r) + 11 
p{r) + 5 



3 r F{r) 



2 \ r 

K.2 



= — 



4 r p(r) + 3 r — r — 1 



(5.20) 



p{r) 



p{r) + 5 
6 r Fir) 



6 r F{r) 



The functions F{r) and p{r) are defined in ( [4. 201) and (|5.17|) , respectively. 

The desired solution to ( ^.19| ) can be found by intergrating the right hand side of the 
equation twice and choosing the solution to the homogenous solution such that we retain 
regularity]^ at r = 1 and appropriate normahzability at infinity. The solution with these 
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We have imposed the requirement that all metric functions are well behaved in its neighbourhood of 



30 



properties is 



dx 



X {x^ 



V. 



dy 2yAij{y) 



(2) 

The quantity of prime interest to us is the leading large r behaviour of a-^. This can 
be inferred from the expressions for the coefficient functions given in ( |5.2(J| ) and evaluates 
to 



r 
1 
2^ 



1 - 



In 2 



(5.21) 



The leading term here will give a divergent contribution to the stress tensor, which is 
necessary to cancel the divergence arising from the expansion of g^^^ + g^^^ to second order. 
The subleading piece in ( p.21| ) is the term that will provide us with the second order stress 
tensor. Before proceeding to evaluate the stress tensor we present the full solution to 
second order, appropriately covariantized. 



5.4 Global solution to second order in derivatives 

Consider the following metric 



ds'^ = -2 dx^dr-r^ f{hr) u^u^ dx^dx^+r^ P^^ dx^dx''+?, h2{br) dx^dr+Q^^ dx^dx" , 

(5.22) 

where we have defined a symmetric tensor ^^j, by combing the contributions in the field 
theory directions from the first and second order metrics f?'-"-* + g'^^^ 



k2{br) 



+ rH^h2{br)P^, + ^ (-2r'Vu^ + ^£\br)] u 



(5.23) 



The covariant expression for ajfj is given by ( |5.19|) with the replacements 



X6,; 



2d, 



{T2,] 



2c, 



2 (T, 



2a, 



3 



(5.24) 



r = 1, a regular point in the spacetime manifold. Note that r = 1 will not represent the horizon of our 
perturbed solution, but may well lie very near this horizon manifold. 
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Further, given by (|5.18| ) with Bj 
make the following replacements 



By{r), where Bj(r) is given by ( ^.15|) and we 



V4i 



7 ^(l3Ua) 



P!f dr.d, 



daU" VUy , 



(5.25) 



232,; 



233. 



Finally, h2{r) and k2{r) are given by ( ^.91) and (|5.11| ) respectively, and in the functions 
Sh{r), Sk{r), S'^ and S'j^ defined in ( ^.8|) and (|5.1CI|) we are required to make the replace- 
ments 



S3 
63 



&2 



64 



©5 



(5.26) 



It may be checked that this metric is the unique (up to terms that differ at third or 
higher order in derivatives) covariant expression that reduces to two derivative solution 
determined in the previous subsections, in the neighbourhood of any point after making 
the coordinate change that sets h^^^ = 1 and /S^^^ = at that point. It follows that ( |5.22|) 
is the desired metric g^^^ + g^^^ + g^'^\ 



5.5 Stress tensor to second order 

The stress tensor dual to the solution to second order described in 5 15.41 can be obtained 



by using the standard formula ( |4.25| ). To determine the extrinsic curvature at large r, it 



suffices to know the asymptotic form of the metric since we are interested in terms that 
have a finite limit as we take r —>■ 00. Consequently, in order to compute the stress tensor it 
is sufficient to replace the various functions of r that have appeared in the computation in 
§ § ^.2| and § ^]3| by their large r asymptotics. The stress tensor may the be computed 
in a straightforward fashion, yielding 

(72)„„ = (72)„. = , 

In 2^ . In2\ / ^ 1^ \ (5-27) 



= -— 17^, - Xe,, + (^-1 + (^t3,, + Ti,, + - 

The vanishing of (^2)^ is actually guaranteed by our renormalization condition. It is easy 



to check that the covariant form of the expression ( ^.27|) is indeed the stress tensor quoted 



in ( p.8|) . This result is the main prediction of our fluctuation analysis to second order in 
the derivative expansion. 
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6 Second order fluid dynamics 



In the previous section we have derived the precise form of the fluid dynamical stress 
tensor dual to gravity on AdSs including all terms with no more than two derivatives. In 
this section we initiate a study of the physics of this stress tensor. In § ^]T| below we will 
demonstrate that our stress tensor transforms homogeneously under Weyl transformations. 



In § we compute the dispersion relation for low frequency sound and shear waves that 



follows from our stress tensor. In § ^]3| we compare our predictions for the transport 



coefficients at second order with those of [37 



6.1 Weyl transformation of the stress tensor 

Thus far we have extracted the stress tensor for a conformal fluid in flat space R^'^. We 
would like to ensure that the second order stress tensor we have derived transforms homo- 
geneously under Weyl rescaling. In order to check this we perform the obvious minimal 
covariantization of our stress tensor to generalize it to a fluid stress tensor about an arbi- 
trary boundary metric g^v^ and study its Weyl transformation properties. 
Consider the Weyl transformation of the boundary metric 



It is well known that the first order truncation of the stress tensor ( p. 81) transforms as 
j^nu _ ^-G^rptiu ^nder this transformation (see for instance Appendix D of [Q)- P^*^" 
ceed to show that this transformation rule holds for the two derivative stress tensor as well. 
This transformation property, together with the tracelessness of the stress tensor, ensures 
Weyl invariance of the fiuid dynamical equations V^T'^'^, appropriate for a conformal fiuid. 

It follows from that P^"' = g^"" + u^u'^ = e'^'f'P^"'. The Christoffel symbols 



transform as 44 



r V = r V + d,<p + 9a0 - gx, r"" 5.0 . 

The transformation of the covariant derivative of u'^ is given by 

= + r;;;, = e-^ [v^ + 6-; u'^ d^<P - g,x -9.0] • (6.2) 
This equation can be used to derive the transformation of various quantities of interest in 



^^AU metrics in this subsection refer to the metric on the boundary, i.e.., the background spacetime on 
which the fluid is propagating. 
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(6.3) 



fluid dynamics, such as the acceleration a'^, shear a'^'^, etc.. 

e = V^u^" = e-^ (Vf.u" + 3 d^(j)^ = e'^ {9 + 3 V(j?j 
a" = Vu" = u^'V^u^ = e-^*^ (a" + P'"' d^<p 

3 

where in the last equation we have accounted for the fact that all epsilon symbols in ( |2.9| ) 
should be generalized in curved space to their covariant counterparts. The objects with 
correct tensor transformation properties scale as metric determinants i.e., eap-ys oc ^Jg, 
and e^^T"^ oc from which it is easy to infer their scaling behaviour under conformal 
transformations; in particular, e^p^s = e^"^ '^aP'^s and e'^^^^ = e~^'^ ^o/375_ 

The Weyl transformation of the two derivative terms that occur in the stress tensor 
]9D is given by 



(6.4) 



Tf = e-^'f' rr , for A = {2a, 26} 

T^- = e-^^ (fj^^-' + Sfs^ , for 5 = {2c, 2c/, 2e} 

where the inhomogeneous terms arising in the Weyl transformation are: 



STJ^" = 2 a^" + 2 u^" a"^ P0 - - a° V« 



2 
3 

3 3 



3 3 

- ^'^(1)0)2 + - P^'^ {V(j)f - 2 n^^ - + - P"^ V„0 

3 3 

(6.5) 

While the conformal transformation involves the inhomogeneous terms presented in 
( |6.5| ) we need to ensure that the full stress tensor is Weyl covariant. Satisfyingly, these 
inhomogeneous terms cancel among themselves in the precise combination that occurs in 
( p.8|) ; consequently the linear combination of terms that occurs in the stress tensor trans- 
forms covariantly. Note that the cancelation of inhomogeneous terms depends sensitively 
on the ratio of coefficients of T2c, and T2e] and so provides a check of our results. Note 
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however that and are separately Weyl covariant. In summary, our result for the 
two derivative stress tensor is a linear combination (with precisely determined coefficients) 
of three independently Weyl covariant forms, with scaling weight —4 (for upper indices). 

Using the transformation of the temperature ( |6.1D it follow that the full stress tensor 
transforms under Weyl transformation as 



rpfll/ 



rpflU 



(6.6) 



R. Loganayagam ^T|] informs us that he has found a compact way of rewriting our stress 



tensor T^'^ (2.S) that makes the Weyl invariance of each of its three pieces manifest. 



6.2 Spectrum of small fluctuations 

Consider a static bath of homogeneous fluid at temperature T. Given the two derivative 
stress tensor derived above ( p.9|) , it is trivial to solve for the spectrum of small oscillations 
of fluid dynamical modes about this background. As the background is translationally 
invariant, these fluctuations can be taken to have the form 

. . , (6.7) 

T{v,x^) = i + STe'^''+'^^'' 

Plugging ( |6.7| ) into the equations of fluid dynamics ( p.7| ), and working to ffist order in 
SPi and 6T, these equations reduce to a set of four homogeneous linear equations in the 
amplitudes and 6T. The coefficients of these equations are functions of u and ki. These 
equations have nontrivial solutions if and only if the matrix formed out of these coefficient 
functions has zero determinant. Setting the determinant of the matrix of coefficients to 
zero one can find the following two dispersion relations: 

Sound mode : oj{k) = ±^ + — ± i^Jli^ + O (k'^) , (6.8) 

v3 6 24 V 3 

• 1 2 

Shear mode : uj{k) = ^ + ^ (2 - In 2) + O (k^) , (6.9) 
where we have defined the rescaled energy and momenta 

71 I 71 I 

It would be interesting to check our prediction against the quasinormal mode analysis 

of m. 
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6.3 Comparison with Baier et.al. (added in v2) 

In this subsection we compare our results with those of the preprint which appeared 
in an arXiv listing simultaneously with the first version of this paper. Where our results 
overlap we find perfect agreement. 

The authors of demonstrated that conformal invariance determines the two deriva- 



tive fluid dynamical stress tensor of any conformal field th eory to upto five parameters 
(see also for an alternate derivation of this result). The five undetermined parameters 



are the coefficients of the various Weyl invariant two derivative expressions built out of 
velocity. These coefficients which are named rn, k and Ai,2,3 by are defined via the 
equation^ 

T,^^^ = fnv rr + ^ rr + Ai + A2 + A3 , (6.11) 



withQ 



2e 



3 3 

= a^^^'a\ = T!^^ , (6.12) 

Here the quantities 6 and a are defined previously (|6.3| ) and oj^'^ = — P^^ <9[a^i/3] is the 
antisymmetric two tensor built from velocity derivatives. In ( |6.12 ) we have also reexpressed 
the operators Tf^ , of [0 as linear combinations of the tensors we have used in 
our paper. These relations are easy to verify given the definitions (p.9|) and (|6.3|). 



The analysis of [B^ was able to determine three of the five coefficients above. Specifi- 



cally, for A/" = 4 Yang Mills they find: 



T-n = n= — , Ai = — , ri=-NT (6.13) 



TT 



The coefficients A2 and A3 were undetermined by their analysis. Translating the results 
of our paper for the second order stress tensor into the language of this subsection and 



^^Note that a''" as defined in (|2J) differs from tliat of Q by a factor of 2 and we define UJ^"' = -f^^". 
We have introduced the tildcd quantities fn etc., to account for the normalization difference. We have 
fn = 2 Tn, Ai = 4 Ai and A2 = —2 A2. 

^^We have used the notation of |^ to indicate the symmetric, transverse tracelessness; for any two 
tensor JF^"^ 

■pit^u) ^ pt^a pu0 j:^^^^ _ 1 pnu pa/3 j:^^ 
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reinstating the overall normalization^ we find 



2 - In 2 

ttT 



2r] ln2 

ttT 



A, = 0. 



(6.14) 



As the coeffient k does not enter the equations of fluid dynamics in flat space, our analysis 
leaves this coefficient undetermined. 

In summary our results for rn and Ai are in agreement with those of Putting 



together the results of our paper with the value of k determined in fS^I we have a pre- 
diction for the full two fluid dynamical stress tensor of a conformal fluid dual to gravity, 
propagating on an arbitrary curved background. 

We can also directly compare the dispersion relations presented in the previous subsec- 
tion with those obtained in 



]. Our sound wave dispersion relation ( |6.8|) agrees with the 
direct computation of quasinormal modes presented in [^. The shear mode dispersion 
relation of the previous subsection also agrees with the direct quasinormal mode analysis 
of upto the order that it should, i.e.. upto terms of order k^. Note that the con- 



tribution to the dispersion of the shear quasinormal mode of [0] does not agree with the 
coefficient of k^ in ( |6.9|) , but (as explained in [^) there is no reason that it should. As 
terms of order k'^ are two orders subleading compared to terms of order k"^ (the order at 
which the first order stress tensor contributes to the shear quazinnormal mode) the coef- 
ficient of /c^ in the shear quasinormal mode potentially receives contributions from third 
order terms in the fluid dynamics stress tensor that we have not accounted for in 



7 Discussion 



We have demonstrated how to start from a general, stationary black brane solution de- 
scribing perfect fluid dynamics and promote the parameters in the gravitational solution to 
physical fluctuation modes. This procedure allows us to set up a fluctuation analysis which 
can be used to extract the boundary stress tensor of fluids dual to gravity in asymptotically 
AdS5 spacetimes, in a derivative expansion. Our procedure is ultralocal: we obtain our 
solution by patching together local tubes of the black brane solution into a global solution 
of Einstein's equations. The fact that our solutions tubewise approximate black branes 
(see 0, 1^, |T^ for related observations) is the gravitational analogue of the fact that the 
fluid dynamics approximation only works when the fluid is in local equilibrium. We find 
this structure of our solutions quite fascinating and feel that it might have the potential 
to teach us important lessons about black brane dynamics. 

^■"^ Recall that in writing ( |2.9| ) we scaled out an overall factor of IGttGn, which evaluates to |^ for 
Af ~ 4 SYM. The value for any other conformal field theory with gravitational dual is simply determined 
by the central charge. 
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Equation ( |2.9| ) is a prediction for the stress tensor of all four dimensional conformal flu- 
ids that admit a dual gravitational description. As we have described in the introduction, 
there exists an infinite number of examples of conformal field theories with a gravitational 
dual that differ substantially in their field content, spectrum of operators, etc. Nonethe- 
less, up to an overall normalization, each of these theories has the same fluid dynamical 
expansion! Consequently, the fluid dynamics described in this paper has a degree of uni- 
versality associated with it. At the one derivative level, the fluid stress tensor has a single 
undetermined parameter - the shear viscosity. The value of rj that we find is in agreement 
with earlier work, rj/s = 1/(4 vr). This relationship has been shown to have a larger degree 
of universality than is apparent from our work; it applies to all field theories, whether 
conformal or not, that have a gravitational dual. This relationship has also been conjec- 
tured to act as a lower bound on the viscosity of a relativistic field theory. It would be 
interesting to investigate whether any of the new two derivative coefficients we have found 
in this paper display extended universality features and also whether they are sensitive 
to higher derivative terms as discussed recently for the shear viscosity to entropy ratio in 



As we have remarked in § 0, it would be interesting to investigate whether our result 
for the stress tensor is consistent with the so called Israel- Stewart formulation of fluid 



dynamics |5T], a framework that has been employed in several practical investigations of 
fluid flows. 



Relatedly, we note that recent claims that the RHIC plasma violate the viscosity 
to entropy bound referred to above are based on the analysis of RHIC plasma flows using 
first order fluid dynamics. However, a satisfactory analysis of these flows should include 
contributions from higher order terms in the fluid dynamical expansion. It is possible that 
the stress tensor derived in this paper will be useful in this regard. ^ 

It may be possible to use the formalism presented in this paper to obtain a better 
understanding of the formal structure of the fluid dynamical expansion of quantum field 
theories. In this context it is useful to recall that the spectrum of regular small oscillations 
about a uniform black brane hosts an infinite spectrum of quasinormal modes. In this 
paper we have effectively constructed the 'chiral Lagrangian' corresponding to those of 
the quasinormal fluctuations that are Goldstone modes (and so have zero frequency when 
at zero k). The remaining quasinormal modes played no role in our analysis, as they 
are nonperturbatively massive in the inverse temperature (tu ~ T = 1/6). The existence 
of these non perturbative modes probably implies that the fluid dynamical expansion is 
asymptotic rather than convergent, and might allow us to predict the location of the first 
singularity in the Borel transform of this perturbation series. 

Recall that metric fluctuation, for any asymptotically AdSs solution to Einstein's equa- 
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We thank O. Aharony for this suggestion. 
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tions, decays at large r like relative to the background. The coefficients of this \/r^ 
decay are functions of the four field theory coordinates in a particular gauge these 
functions may be identified with the 9 components of the traceless boundary stress tensor. 
This stress tensor is constrained to obey the equations of energy momentum conserva- 
tion, but is otherwise unconstrained by local analysis. The Fefferman-Graham |^ method 
(or equivalently the formalism of holographic renormalization, see |5D[ for reviews) 
demonstrates that any such conserved stress tensor, regarded as a boundary condition to 
Einstein's equations, leads to a unique and well defined power series expansion (in 1/r) of 
an asymptotically AdS metric. Local analysis near the boundary thus appears to indicate 
that the space of solutions to Einstein's equations in AdS space is parameterized by the 
set of all conserved energy momentum tensors in four dimensions. This would be very 
surprising from the dual field theory viewpoint, as a set of four equations does not define 
a well posed initial value problem for nine functions. 

The results of our paper suggest a (perhaps not unanticipated) resolution to this puzzle. 
In the derivative expansion in which we work, all except a four function set of this naive 
nine function class of metrics are unacceptably singular and so do not constitute a legal 
solution to Einstein's equations. Generic data result in singularities that develop at a finite 
value of r (r = 1/6 in our set up) and so are not easily visible in the Fefferman-Graham 
expansion, which is guaranteed to work only in an open neighbourhood of the boundary. 
The class of boundary stress tensors that generate acceptable metrics are parameterized 
by four functions {(3i{x^) and h{x^)) rather than nine. These four functions are further 
constrained to obey the four equations of stress energy conservation. As four equations 
constitute a well defined^ initial value problem for a set of four functions, the set of 
legal solutions to Einstein's equations are parameterized by data that consists of functions 
of 3 spatial rather than 4 spacetime boundary variables, in agreement with field theory 
expectations. It would of course be of very great interest to understand how these results 
of the previous paragraph generalize beyond the boundary derivative expansion. 

In this context it is also relevant to note that the equations of fluid dynamics themselves 
develop singularities under certain situations. It would be interesting to investigate the 
gravitational dual of this process of singularity formation.^ More generally, the map 
from solutions of fluid dynamics to solutions of gravity could allow one to use the insight 
gained from the hundred year long study of the equations of fluid dynamics to understand 
qualitatively new gravitational solutions. For example, one might hope to learn about 
stationary inhomogeneous brane solutions (analogous to those discovered in the study 



^''Note that our notion of a well posed PDE system is simply that we do not have an under-constrained 
system of equations. We are not making any claim regarding the well posedness of generic initial data; only 
initial data in the regime of our perturbation analysis together with the boundary conditions is guaranteed 
to lead to regular solutions. The general question of global regularity of Navicr-Stokcs equation is of course 
an interesting open problem. 

^®We thank D. Berenstein for discussions on this question. 
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of Gregory-Laflamme instabilities of black strings and branes) using the fluid dynamical 
description. 

As we have described, in this paper we have derived explicit formulae for the metric dual 
to any solution of the Navier-Stokes equations. We have not yet investigated the global 
structure of the resulting spacetime. It seems very plausible that (under suitable physical 
conditions) the spacetimes we have constructed have regular event horizons. The event 
horizon is a null surface; we expect it to closely approximate^ the surface rh{x^) = 1. 
If this is the case we should be able to compute an explicit expression for this surface 
order by order in perturbation theory. It may then be possible to use our understanding 
of the horizon to define a locally positive divergence entropy current (a 'pullback' of the 
one-form dual to the natural null generators of the horizon back to the boundary might 
play a role in such a construction). In the most optimistic scenario such an exercise could 
relate classic results about the positivity of null congruence expansions (resulting from the 
Raychaudhuri's equation with the usual proviso of energy conditions) to the local positivity 
of entropy production in fluid dynamics; a result that would be of obvious interest. The 
language of dynamical horizons^ may well prove to be the appropriate framework for 
such a discussion. We hope to return to this intriguing issue in the future. 

Recall that the construction presented in this paper yields the gravitational dual of 
every solution of the equations of fluid dynamics. Standard field theory lore asserts that 
generic field theory evolutions are well described by solutions to the equations of fluid 
dynamics in the regime of interest to this paper. Consequently the AdS/CFT correspon- 
dence seems to imply that the construction described here yields the generic legal solution 
of gravity in AdSs, within its domain of applicability. If our guess of the previous para- 
graph is correct - i.e., if all our solutions possess a regular event horizon that shield the 
boundary from the singularity - then our results appear to be of relevance to the cosmic 
censorship conjecture p2| , |53| . However we emphasize that our analysis applies only in a 
long wavelength expansion, and so presumably does not apply to several classes of scenarios 
that putatively violate this conjecture. More physically, fluid dynamics applies only under 
the assumption of local thermal equilibrium. Presumably naked singularities (if they exist) 
are dual to 'far from local equilibrium' boundary physics. 

Several other natural generalizations of the work presented in this paper immediately 
suggest themselves. First, the results of our paper are likely to have an analogue for d 
dimensional gravitational theories with a negative cosmological constant for every d > A. 
Second, it should be possible to extend the results of our paper to spaces that asymptote 
to AdSd+i in global coordinates (and whose dual description is, therefore, fluid dynamics 

^^We would like to thank H. Reall for many useful discussions on this point. 

■^°Note that the homogeneous spacetime background being simply the uniform black brane, has a trapped 
surface; under the fluctuations we generically expect the trapped surface to be generated earlier in the 
radial evolution i.e., at a larger radial coordinate (assuming of course appropriate energy conditions). 
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on S'^ ^ X Rf). More ambitiously, it may be possible to extend the results of our work to 
field theories whose spacetime metric asymptotes to 

for a more general class of metrics dsf^y. In particular, the generalization to time dependent 
metrics would permit the study of the gravitational dual of forced fluids, a subject of 
interest to the study of turbulence. Finally, it should not be difficult to generalize our 
study to two derivative theories of gravity interacting with gauge fields. We expect that 
the dual description of such a system will be the fluid dynamics of a system with a number 
of additional conserved charges (equal to the number of commuting vector fields). Note 
however that unlike the uncharged system, this charged fluid dynamics will not be universal 
at nonlinear order, as gauged supergravities do not in general admit a consistent truncation 
to the Einstein-Maxwell sector. For instance couplings of the form f{(j))F^^F^^ , for an 
arbitrary scalar field 0, constitute a source for </>; this is an effect that plays an important 
role in studies of the attractor mechanism. 

Despite this non universality, IIB supergravity on AdSs xS*^ (for instance) should be 
dual to a completely well defined charged fiuid dynamics. It would be of interest to use 
the methods of this paper determine the form of this fiuid dynamical stress tensor. Among 
other things, this exercise would allow us to zero in on the origin of the worrying apparent 
discrepancy between the formulas of charged black hole thermodynamics and the formulas 
of fluid dynamics, as reported in f^. 

It should prove relatively straightforward to generalize the study of this paper to the 
fluid dynamics of non conformal backgrounds of gravity. For example, Scherk-Schwarz 
compactiflcations of AdS spaces yield a particularly simple set of gravitational backgrounds 
dual to confining gauge theories. Indeed, a moment's thought is enough to convince oneself 
that the deconfined phase fiuid dynamical stress tensor of the 2+1 dimensional confining 
gauge theory (dual to Scherk-Schwarz compactified A/" = 4 Yang Mills) is simply the 
dimensional reduction of the stress tensor of d = 4, A/" = 4 Yang Mills {i.e., the stress tensor 
presented in this paper, ( |2.8| )) plus a constant additive piece. This seemingly trival additive 
piece is physically very important; it leads to qualitatively new phenomena. For instance, 
5^ , |55i have plasmaballs and plasmarings; static finite lumps of fiuid with a boundary. Such 



configurations have qualitatively new classes of excitations; localized collective coordinates 
associated with fiuctuations of the boundary. These new collective coordinates will interact 
with those studied in this paper. If it proves to be technically possible, it would be 
fascinating to formulate and study the resulting dynamics. 
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